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Abstract 

The static vacuum spherically symmetric solutions of massive gravity theories possess two inte- 
gration constant: the mass M and a scalar charge S. The presence of this scalar charge reflects 
the modification of the gravitational interaction as compared to General Relativity. Surprisingly, 
these solutions have an asymptotic behavior different from that obtained in the linear perturbation 
theory. The aim of this paper is to understand how these modified spherically symmetric solutions 
emerge from a quasi-linear approximation in order to generalize them to any arbitrary mass dis- 
tribution. Along with these modified solutions, we found a new class of static solutions having a 
Yukawa shape. 
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I — ^1 Introduction 

o 

I Gravity is the most familiar but also the most mysterious of 

H„the four fundamental interactions. Described by Einstein's 

i__i theory of General Relativity (GR), gravity has been tested 

in the weak field limit up to post- Newtonian corrections [1]. 

^^ Although indirect detection of gravitational waves seems to 

J^ be in agreement with GR [2, 3], recent advances in observa- 

^^■ tional cosmology [4, 5, 6] have revived interest in alternative 

l/^ theories of gravity in which the gravitational interaction is 

■^ modified at very large distances. The motivations behind 

l/~j these theories lie in the possibility of explaining the observa- 

^^ tions without introducing the dark energy and matter com- 

^D ponents. Different approaches to this problem have been 

^. discussed in the literature [7, 8, 9, 10, 11, 12, 13, 14]. One 

t> of them employs spontaneous breaking of Lorentz symmetry 

by space-time dependent condensates of scalar fields [15, 16] 

. : in order to give a mass to gravitons^, which is the reason to 

C^ ball these models massive gravity models. 

If the gravitons were to have a mass, one would expect the 
potential of a static source to have a Yuliawa shape effec- 
tively cutting off the gravitational interaction at distances 
larger than the inverse graviton mass ]19, 20, 21]. Absence of 
Lorentz invariance makes these constraints weaker than one 
would expect in a Lorentz-invariant theory; there is a class 
of massive gravity models in which the Newtonian potential 
remains unmodified in the linear approximation ]22], im- 
plying that Solar system constraints are satisfied for rather 
large graviton mass. 

Recently, a class of static vacuum spherically symmetric 
solutions in massive gravity theories were obtained both an- 
alytically and numerically ]23]. These solutions depend on 
two integration constants instead of one in GR: the mass 
M (or, equivalently, the Schwarzschild radius) and an ad- 
ditional scalar charge S. Depending on the parameters of 
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^See Refs. [17, 18] for reviews. 



the model considered and the signs of Af and S, these mod- 
ified solutions may have both attractive and repulsive (anti- 
gravitating) behavior at large distances. Those with at- 
tractive behavior may mimics the presence of dark matter. 
Another peculiarity of these solutions is that their asymp- 
totic behavior is different from that obtained in the linear 
perturbation theory. Despite the obvious interest of these 
solutions, the underlying hypothesis of spherical symmetry 
makes their use limited. 

The aim of this paper is to understand how these modi- 
fied solutions emerge from a quasi-linear approximation of 
massive gravity theories. Beside pure curiosity, the answer 
to this question sheds lights on the shape of the Newto- 
nian potential produced by any arbitrary mass distribution 
by identifying the mechanism and equations responsible for 
these modified solutions. One will then be able to determine, 
for example, the gravitational potential in a galaxy which 
is made of two contributions: one contribution coming from 
the spherically symmetric bulge, and one coming from the 
axially symmetric galaxy disk. Because of the non-linearity 
of the equations responsible for these modified solutions, 
these two contributions cannot be considered separately as 
in linearized GR. 

By studying the perturbations up to second order about 
the flat vacuum solution, we found two classes of static 
spherically symmetric solutions. The first one is a new class 
of solutions corresponding to solutions of a four-derivative 
equation for the Newtonian potential. They have a Yukawa 
shape effectively cutting off the gravitational interaction at 
distances larger than the inverse graviton mass. The second 
class of solutions correspond to the modified solutions found 
in Ref. ]23]. 

This paper is organized as follows. In Sect. 2 we briefly 
review the massive gravity models and summarize previous 
results about the Newtonian approximation and the spheri- 
cally symmetric vacuum solutions. In Sect. 3, we introduce 
notations and discuss the quasi-linear equations which will 



enable us to recover the modified spherically symmetric vac- 
uum solutions. Finally, in Sect. 4 we discuss the two classes 
of solutions. For the second class, we show that our quasi- 
linear static equations reproduce the modified Newtonian 
potential obtained in Ref. [23]. Sect. 5 contains the sum- 
mary and discussion of our results. 

2 The massive gravity model 

In this paper, we consider massive gravity models described 
by the following action [16], 
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Crn+A^J'{X,W'^)] 



(1) 



The first two terms are the usual Einstein-Hilbert term and 
the Lagrangian of the minimally coupled ordinary matter; 
they comprise the standard GR action. The third term de- 
scribes four scalar fields (f)^ whose mixing with the met- 
ric give a mass to gravitons. These fields, known as the 
Goldstone fields, are minimally coupled to gravity through 
a derivative coupling 
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g^^-'d^c^^d.^cb", 



flV 



X 



dn(b'd^CJ)'. 



The Goldstone fields (p^ have dimension of length so that X 
and W^^ are dimensionless. The constant A is a UV cutoff 
with dimension of mass. Massive gravity models (1) are 
understood as low-energy effective theories valid for energies 
below A. 

The Goldstone fields have space-time dependent vacuum 
expectation values that break spontaneously the Lorentz 
symmetry 



gfj^i' = Vfii^^ 4>° = at, cj)" = 6a;' 



(2) 



The constants a and b are determined by the requirement 
that the energy-momentum tensor associated with the four 
scalar fields vanishes in Minkowski background. Once a 
function F is chosen, these constants may be set to one 
by a redefinition of the fields 0'^ . which we assume to be the 
case in what follows, a = b = 1. 

We consider only functions T which are invariant under 
rotations of the Goldstone fields (jf in the internal space 
(i.e., functions depending on W''^ through three combina- 
tions Wn = Tr {W^) with n = 1,2, 3). For these functions, 
the background (2) preserves the rotational symmetry. 

Finally, the Goldstone action with functions T depending 
only on X and W^^ is invariant under the following symme- 
try 

where x* are arbitrary functions of (/)°. Because of this sym- 
metry, the perturbations about the vacuum solution (2) are 
non-pathological, i.e. there are neither ghost nor rapid in- 
stabilities ]16]. The spectrum consist of two tensor modes 
(graviton polarizations) only, which are, in general, massive. 
The graviton mass scale is tti ~ A^/Mpi. 



Before discussing quasi-linear static solutions, let us 
briefiy resume the previous results concerning the Newto- 
nian approximation ]22] and the Schwarzschild solutions [23] 
of theories described by the action (1). 

Newtonian approximation 

For massive gravity models described by (1), the gravita- 
tional potential of a static source, in the linearized approxi- 
mation, has the form [22] 



^ = MG ( — +/^^r 
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where G 
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is Newton's constant. The first con- 



tribution to this potential is the usual GR term. The second 
contribution is proportional to a constant /i which is of order 
of the graviton mass and whose value depends on the partic- 
ular form of the function T. This contribution implies that 
the linearized theory breaks at distances r ^ {MGfjL^^ 

Spherically symmetric solutions 

Recently, the static vacuum spherically symmetric solutions 
of massive gravity theories described by (1) have been ob- 
tained both analytically and numerically [23] for different 
functions J- . In order to get some insight into the behavior 
of these exact solutions, Ref. [23] studies the exact solution 
of models described by the following function 
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where A > is an arbitrary positive constant. For this 
class of models, the exact static solution of the equations of 
motion implies that Newton's potential is 



q> ^-- 



GM 
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where S is an integration constant whose presence reflects 
the modiflcation of the gravitational interaction as com- 
pared to GR. Similar solutions have been previously found 
in the context of bi-gravity models [24]. This potential has 
an asymptotic behavior different from that obtained in the 
linear perturbation theory. The origin of this apparent con- 
tradiction is twofold. 

In the linear perturbation theory, all perturbations are as- 
sumed to be of the same order: they are formally assigned a 
small parameter e to the first power and the Einstein equa- 
tions are expanded to the linear order in e. The modified 
solution (5) is not of this type, even at large distances from 
the gravitational source. For instance, in the longitudinal 
gauge, the perturbations of the Goldstone field (fP does not 
decay as fast as the perturbations of the metric; if one as- 
signs a small parameter e to the perturbation of goo, then 
the perturbation of 0" is of order y/e rather than e. There- 
fore, the exact solution of the static, spherically symmetric, 
massive gravitational field equations is said to be non-linear. 



Another difference between the exact solution (5) and the 
solution (3) of the linearized equations is that the former is 
static, while in the latter only gauge-invariant metric per- 
turbations are static [22]; the perturbation of the Goldstone 
field (jp only appears in the linearized theory with a time 
derivative which may be viewed as an accretion of a fluid 
with zero energy-momentum tensor. Therefore, the two so- 
lutions (3) and (5) simply describe different physical con- 
texts. 

3 Quasi-linear approximation 

Despite the obvious interest of the modified solution (5), the 
underlying hypothesis of spherical symmetry makes its use 
limited; it cannot be used to describe the potential induced 
by an axially symmetric galaxy. It is therefore important 
to understand how this modified solution emerges from a 
quasi-linear theory in order to generalize it to any mass dis- 
tribution. 

To achieve this, we consider perturbations of the metric 
and Goldstone fields up to second order about the flat vac- 
uum solution (2) 

2 2 



and four scalar gauge-invariant fields 
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where |/i(a)^ii/| ^ \^u)\ ^ ^''' '^ith a = 1,2. These perturba- 
tions are parameterized as the following [25] 

h(a)QO = "^fia), 






(6) 



The vector perturbations Sla)i^ ^(a)i s-nd CLi ^^^ trans- 
verse, while the tensor perturbations Hu\^j are transverse 
and traceless. In order to simplify the discussion of second 
order perturbations, let us fix the gauge at the linearized 
level to be the longitudinal gauge 



B{i) — Ell) — F,i 



(1)' 
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This choice for the linearized approximation will greatly sim- 
plify the coming calculations without fixing completely the 
gauge in the perturbation theory. One is still allowed to 
perform second order gauge transformations 



x"^ -> x^ 



r 



with 
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which leave the theory unchanged ]26] . It will then be useful 
to introduce second order gauge-invariant fields. One vector 
and two scalar perturbations are gauge degrees of freedom. 
As a consequence, at the quadratic order there is only two 
gauge-invariant vector fields 
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2° = ^("a) - -8(2) - £^(2) , S - C(2) - E^2) ■ 

The tensor perturbation H(2)ij is also gauge invariant. In 
order to increase the readability of the coming calculations, 
the subscript (a) will be omitted from this point without 
possible confusion between first and second order perturba- 
tions. 

Usual matter fields are described by an energy-momentum 
tensor T^^v It is convenient to parameterize this tensor in 
the following way 



' iif 



{p + p) v^v^ - g^^p + TT^^ + (v^q^ + v^q^) , 
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where p and p are the matter density and pressure mea- 
sured by a comoving observer of velocity u^, g^ is the en- 
ergy fiux perpendicular to w^ {v^^q^ = 0) and tt^^ is the 
anisotropic pressure tensor {v^tt^v = Tri^ = 0). For the vac- 
uum solution (2), the energy- momentum tensor is zero and 
the affine parameter of the observer may be chosen such that 
z;^ = (1,0,0,0). 

Small perturbations 6Tfj,i, of the matter energy - momen- 
tum tensor produce metric (5.9^1/ — g^u ~ rj^u and Gold- 
stone 5(j)^ = (pf^ — x^ perturbations about the vacuum. In 
other words, the metric and Goldstone perturbations are 
sourced by 5T^^. Since the modified solution (5) behaves 
as If^^^p ^ \5gm\ ~ M/r, let us assume that the perturba- 
tions of the gravitational source appear at second order only, 
i.e., \6p\ ^ \5p\ ~ \5qfj\ ^ \5iT^^\ ^ e^. As a consequence, 
(5(7o = (^TToi/ = and the second order energy - momentum 
tensor reads 

<57oo = Sp, SToi = 5qt, STij = 5,,jSp + JTr^ . 

The energy flux Sqi and the anisotropic pressure Sttij can be 
parameterized in the following way, 

6qi = Ci + diC, 
Sn^j = (3didj - Sijdl) n + d.,Tij + OjTt^ + n^j , 

where the vector perturbations Q and tt^ are transverse, 
while the tensor perturbation TTy is transverse and trace- 
less. 

With these notations, the energy-momentum conservation 
reads 

Sp^dfC, C^5p + 2df7T, Q=d]n,. 

The ten independent components of the energy-momentum 
tensor are expressed trough four gauge-invariant scalars 5p, 
Sp, C and TT, four vector degrees of freedom in the form of 
two transverse gauge-invariant vectors Q and tt^, and two 
tensor degrees of freedom represented by the transverse and 
traceless gauge- invariant tensor tt^. 

3.1 First order perturbations 

With the account of the previous notations for the metric 
and Goldstone perturbations, the linearized Einstein equa- 



tioiis split into scalar, vector and tensor equations. These 
tree sectors could be considered separately^. 

Without matter fields at the linearized level to source the 
tensor and vector perturbations, the solutions to the first 
order equations in these two sectors reads 



e 



H, 



0. 



(7) 



The key observation which will enable us to recover the 
exact, static, potential (5) from a quasi-linear approxima- 
tion lies in the scalar sector. Since the perturbation ^" of 
the Goldstone field (fp only appear at the linearized level 
with a time derivative, the solution to the linearized scalar 
equations reads (recall that there is no first order matter 
source) 



The vector sector 

The second order vector perturbations satisfy the two fol- 
lowing relations 



= dlwi — 
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(9) 



The first of these two equations gives the metric perturba- 
tion zui while the second fixes Ci once ^'^ is known. 



The scalar sector 



¥^ = V = e = e" = 0, 



(8) 



implying that ^° is an arbitrary function of the space coor- 
dinates ^° = ^"^ i^'')- This function is not fixed by the first 
order equations, although one usually assumes from them 
that ^^ = 0. Still, it may be that this arbitrary function is 
determined by higher-order equations. Let us therefore not 
fix it at this stage, and study the second order approxima- 
tion. 



3.2 Second order perturbations 

In order to discuss the second order equations, let us as- 
sume that the first order solution is given by (7) and (8) 
with ^^ (a:*) to be determined. As for the linearized equa- 
tions, the second order equations split into scalar, vector 
and tensor equations. However, each of these three sectors 
contain contributions proportional to ^q. 



The second order scalar sector is described by the following 
four equations 



= 
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(10) 

(11) 

[3* + dj^ 
(12) 



(13) 



The tensor sector 

The second order tensor perturbations satisfy the following 
equation 



(9^5^ + ml) H,, 



2tt, 
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where 1712 c>c A^/Mpi is the graviton mass (there are five 
mass parameters rrii with i = 0,...,4; see Ref. [27] for a 
precise expression of these masses in terms of T{X, W^^ ) 
and its derivatives). This equation describes gravitational 
waves Hij emitted by the Goldstone perturbation ^° and 
a source with transverse and traceless anisotropic pressure 

TTij. 



^See the appendix of Ref. [27] for details about the linearized equa- 
tions of theories described by (1). 



This system of four equations describe the four potential $, 
*, S and S° emitted by the Goldstone field $° = $° (2:*) and 
a source with matter and pressure densities Sp and Sp. 

At this point, the ten second order Einstein equations 
have been written and one has to conclude that ^° is not 
determined by second order equations. Moreover, the sec- 
ond order equations do only constraint the value of 'EP , and 
therefore the value of S° modulo an arbitrary function of 
the space coordinates. Consequently, the third-order equa- 
tions contain ten independent, third-order, gauge-invariant 
fields'^ plus two arbitrary functions of the space coordinates 
corresponding to the first and second order initial conditions 
on (fP. Hence, these two arbitrary functions wont be con- 
strained by third-order equations either. One may convince 
oneself that, at each order, the Einstein equations determine 
0° modulo an arbitrary function which is not fixed by any 
equation. This arbitrary function may then be chosen to be 



^There are six third-order, gauge-invariant, fields coming from the 
"pure" gravitational sector and four from the Goldstone sector. 

*The addition of higher-derivative terms in the theory should lift 
the ambiguity concerning these arbitrary functions. 



4 Static quasi-linear solutions 

Although the previous argument seems valid for any arbi- 
trary field configuration, it has a loophole. Indeed, let us 
search for static solutions to the second order equations, 
with 



whose solution for a point-like source 5p = M5'^ {x) may be 
approximated at small and large distances by 

MG 
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for 
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Sp — Sp (a;*) and 6p — 0. 

Then, equation (11) reduces to an identity and S" disap- 
pears from the equations of motions''. One is therefore left 
with three scalar second order equations for ^, $ and S, 
and the second order Einstein equations may be solved by 
assuming ^'' = 0. However, by doing this one will miss the 
branch of solutions which reproduces the exact, modified 
Newtonian potential (5). 

For static solutions, there is actually a fourth equation 
which as to be satisfied by second order perturbations. This 
equation correspond to the longitudinal part of the third- 
order approximation of the following Einstein equation 
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where Q^^i^ is the Einstein tensor and t^^, is the energy- 
momentum tensor of the four Goldstone fields. Indeed, for 
a static solution, this third-order equation reduces to 
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Since there is no third-order fields in this equation^', any so- 
lution to the second order equations has to satisfy this third- 
order equation as well. Moreover, this equation implies that 
there are two branches of static solutions, as already stressed 
in Ref. [23]. The first branch corresponds to solutions with 
^° = for which equation (14) reduces to an identity. The 
second branch, characterized by ^° 7^ 0, is the one studied 
in Ref. [23]. Let us discuss these two branches separately. 



4.1 Solutions with ^^■ 







If one considers solutions for which ^° = 0, then the second 
order scalar equations reduce to a system of three equations 
(10), (12) and (13) for *, $ and S. These equations imply 
that Newton's potential satisfies a four-derivative equation 
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^This field only appear in the second order Einstein equation with 
a time-derivative. 

®The third-order fields enter this equation with a time-derivative 
and therefore disappear when considering static solutions. 



where rj and r are dimensionless combinations of the masses 
777-1 . One therefore concludes that the Newtonian potential 
deduced from a static solutions with ^° = is given by its 
value in GR at small distances and is cut off at distances 
larger than 777"^. 

4.2 Solutions with ^° ^ 

The solutions with ^" 7^ are those which reproduce the 
modified potential (5). Because of the complexity of the 
third-order equation (14), let us concentrate on spherically 
symmetric solutions. 

For static, spherically symmetric solutions with ^^ 7^ 0, 
equation (14) reduces to 



= Uml-ml-ml + 2m^:){dref 
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while the three equations (10), (12) and (13) read 
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This is a system of four equations for $, ^I*, S and (5,-^o) • 
It can be solved as follows. Adding equations (16), (17) and 
(18) to equation (15) gives a closed equation for $ + rdr^ 



= ^-dr ($ + 7-a,*) - ^ 



(19) 



(20) 



whose solution reads 
$ + 79^* = w(r), 
w (r) = 47rG 



For a point-like source 8p — M6^ (x) and w (r) = —GM5 (r), 
where 6 (r) is the delta function in spherical coordinates ' . It 



dr' / dr rS (r — r') Sp [r') 



'^The one dimensional delta function 5 (r) is related to the three 
dimensional one through the usual relation 
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is worth noting that equation (19) is the linearized version of 
the relation dr{goo grr) = characterizing the Schwarzschild 
solution of GR in the usual spherical coordinates. 

Solving the three remaining equations for arbitrary func- 
tions J^ of X and W^^ is a difficult task which will not be 
achieved here. However, let us discuss the solution of these 
three equations for two classes of models. 

First class of models 

The first class of models which will be addressed here corre- 
spond to models described by the function (4) discussed in 
Ref. [23] . We should stress here that this class of models is 
only particular because we know how to solve analytically 
the exact, static, spherically symmetric equations in these 
models. 

At the quasi-linear level, the static equations are easily 
solved because of the following relations between the mass 
parameters 

niQ = 3ml, 1712 = "*4 (2 + A) , rn^ = ml (1 + A) . 

Indeed, having the solution of equation (19) in nund, equa- 
tions (15), (16) and (17) imply that 

dr - ^ {drE + r*) = w (r) . 



The solution to this equation reads 
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For a point-like source, k (r) = 0. The remaining two equa- 
tions imply then a closed inhomogeneous equation for the 
Newtonian potential $ 
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For a point-like source, the solution of this equation is ex- 
actly the modified potential (5) 
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where r^ and S are two integration constants. Making use 
of this solution and integrating equations (20), (21) and (18) 
one find the three other fields 
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Second class of models 

The second class of models which will be addressed here 
correspond to attractors of the cosmological evolution [22]. 
Their action depends on the variables X and W^^ through 
a single variable Z*-' = X'^W'^^ 



F = F(Z'^] 



(23) 



where the constant 7 is a free parameter. These models 
have been studied intensively [22, 27, 28]. In particular, it 
has been shown that for —1 < 7 < and for 7 = 1 the cos- 
mological perturbations in these models behave identically 
to those in GR [27]. For other values of 7 the behavior of 
the cosmological perturbations may or may not reproduce 
that of GR depending on the initial conditions. 

Another reason to study models characterized by a func- 
tion (23) comes from the analysis of the linearized approxi- 
mation. For these models, the /z parameter in equation (3) is 
zero [22[ with consequences that Newton's potential remains 
unmodified at the linearized level. 

Still, at the non-linear level, the spherically symmetric 
vacuum solutions have similar behavior [23] as those de- 
scribed by the modified potential (5). Indeed, equations 
(15), (16) and (17) are easily solved because of the following 
relations between the mass parameters 



TOn 



37771; 



7 (3m; 



characterizing these models. Having the solution of equation 
(19) in mind, equations (15), (16) and (17) imply that the 
Newtonian potential satisfies 
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Here, A = 1 — I/7. For a point-like source, the solution of 
this equation is also given by the modified potential (22). 
The behavior of this solution at r -^ 00 depends on the 
constant 7. For 7 < and for 7 > 1, this solution describes 
asymptotically flat space-time. 

5 Discussion 

The aim of this paper was to understand how the mod- 
ified Newtonian potential (5), inferred from static spheri- 
cally symmetric solutions of massive gravity theories, arises 
from a quasi-linear perturbation theory. By studying the 
Einstein equations up to second order perturbations, one 
has identified two branches of static solutions. In the first- 
order longitudinal gauge, the key difference between these 
two branches is the perturbation ^" of the Goldstone scalar 
field 0°. 

The first branch of solutions is a new one corresponding to 
solutions with ^° = 0. They are solutions of a four-derivative 
equation for the Newtonian potential. They reproduce the 
potential of GR at small distances with a Yukawa shape ef- 
fectively cutting off the gravitational interaction at distances 
larger than m~^ . 

The second branch of solutions, characterized by ^^ 7^ 0, 
are those who reproduce the modified potential (5). The 



quasi-linear equations have been solved with ^'^ 7^ in two 
different classes of models. The first one corresponds to 
models characterized by the function (4). The static vac- 
uum spherically symmetric solutions were already known in 
these models. The second class of models correspond to at- 
tractor of the cosmological evolution for which one only has 
numerical static spherically symmetric solutions. By solv- 
ing the quasi-linear equations, the present work sheds lights 
on the static spherically symmetric solutions in this second 
class of models. 

In this second branch of solutions, the Newtonian po- 
tential (22) depends on two integration constant: the 
Schwarzschild radius r^ and an additional scalar charge S. 
For usual source, these two integration constants have to be 
zero. Indeed, for a sphere of constant density g and radius 
R, the energy density is given by 6p = qQ {R ~ r). For such 
source, one may show that the solution to equation (24) 
reads 



$ = 



3GM 
2R 

GM 



GM r, S 
Q{r- R) J 



2R 



r 

i?2 



3 e (i? - r) 



with M = AttqR^ /?,. The first term in this equation is a 
new contribution which may be re-absorbed by a redefinition 
of time. Since there is no singularity in the center of the 
source, the only physical solution corresponds to r^ = 5 = 0. 
Therefore, the scalar charge of an ordinary source is zero. 

It remains an open question how objects with S ^ Q can 
be created. However, the argument given above does not 
apply to black holes, especially to super-massive black holes 
in the centers of galaxies, which may be of primordial ori- 
gin [29]. The previous argument does not apply to time- 
dependent configurations either, so it is possible that a non- 
zero scalar charge may be acquired during the gravitational 
collapse. 

To conclude, let us stress that the equations presented in 
this paper may be used to determine the static gravitational 
potential produced by any arbitrary mass distribution. Al- 
though we have mainly discussed spherically symmetric so- 
lutions, the potential produced by a galaxy, for example, 
may be found by solving numerically equations (9), (10), 
(12), (13) and (14) altogether with the hypothesis of axial 
symmetry. 
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